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Steady scale-invariant solutions of a kinetic equation describing the statistics of oceanic internal 
gravity waves based on wave turbulence theory are investigated. It is shown that the scale-invariant 
collision integral in the kinetic equation diverges for almost all spectral power-law exponents. These 
divergences come from resonant interactions with infra-red and/or ultra-violet wavenumbers with 
extreme scale-separations. 

We identify a small domain in which the scale-invariant collision integral converges and numer- 
ically find a convergent power-law solution. This numerical solution is close to the Garrett-Munk 
spectrum. Power-law exponents which potentially permit a balance between the infra-red and ultra- 
violet divergences are also discussed. The balanced exponents are generalizations of an exact solution 
of the scale-invariant kinetic equation, the Pelinovsky-Raevsky spectrum. 

Ways to regularize this divergence are considered. A small but finite Coriolis parameter repre- 
senting the effects of rotation is introduced into the kinetic equation to determine solutions over the 
divergent part of the domain using rigorous asymptotic arguments. This gives rise to the induced 
diffusion regime. 

The derivation of the kinetic equation is based on an assumption of weak nonlinearity. Dominance 
of the nonlocal interactions puts the self-consistency of the kinetic equation at risk. Yet these weakly 
nonlinear stationary states are consistent with much of the observational evidence. 

PACS numbers: 47.35.Bb 



I. INTRODUCTION 



Wave-wave interactions in continuously stratified flows have been a subject of intensive research in the last few 
decades. Of particular importance is the observation of a nearly universal internal-wave energy spectrum in the 
ocean, first described by Garrett and Munk. It is generally thought that the existence of a universal spectrum is 
at least partly the result of nonlinear interactions among internal waves. Due to the quadratic nonlinearity of the 
underlying fluid equations and dispersion relation allowing three-wave resonances, internal waves interact through 
triads. In the weakly nonlinear regime, the nonlinear interactions among internal waves concentrate on their resonant 
set, and can be described by a kinetic equation, which assumes the familiar fonnii2AiiSi£ffl i 10 i 1 1 i 12 i 13 t 14 : 



_\ypi |2 fpi s: r 

\ v P2,p\ JP2P U P1~P2-P "u P1 -a) P2 -iii, 
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with fP ip2 = n Pl n P2 - n p (n Pl + n P2 ) . (1) 

Here n p = n(p) is a three-dimensional action spectrum (see Eq. (|16p ) with wavenumber p — (k,m), i.e. k and m 
are the horizontal and vertical components of p. The frequency lo is given by a linear dispersion relation. The 
wavenumbers are signed variables, while the wave frequencies are always positive. The factor V pip2 is the interaction 
matrix describing the transfer of wave action among the members of a resonant triad. 

Following Kolmogorov's viewpoint of energy cascades in isotropic Navier-Stokes turbulence, one may look for 
statistically stationary states using scale-invariant solutions to the kinetic equation JTJ) . The solution may occur in an 
inertial subrange of wavenumbers and frequencies that are far from those where forcing and dissipation act, and also 
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far from characteristic scales of the system, including the Coriolis frequency resulting from the rotation of the Earth, 
the buoyancy frequency due to stratification and the ocean depth. Under these assumptions, the dispersion relation 
and the interaction matrix elements are locally scale- invariant. It is natural, therefore, in this restricted domain, to 
look for self-similar solutions of Eq. (UJ) , which take the form 

n{k,m) = \k\- a \m\- b . (2) 

Values of a and b in the right-hand side of Eq. ([1]) vanishes identically correspond to steady solutions of the kinetic 
equation, and hopefully also to statistically steady states of the ocean's wave field. Unlike Kolmogorov turbulence, 
the exponents which give steady solutions can not be determined by the dimensional analysis alone (see, for example, 
Ref. [la). This is the case owing to multiple characteristic length scales in anisotropic systems. 

Before seeking steady solutions, however, one should find out whether the improper integrals in the kinetic equa- 
tion (fT]) converge. This is related to the question of locality of the interactions: a convergent integral characterizes 
the physical scenario where interactions of neighboring wavenumbers dominate the evolution of the wave spectrum, 
while a divergent one implies that distant, nonlocal interactions in the wavenumber space dominate. 

It turns out that the internal-wave collision integral diverges for almost all values of a and b. In particular, the 
collision integral has an infra-red (IR) divergence at zero, i.e. |fci| or |fc2| — * and an ultra-violet (UV) divergence at 
infinity, i.e. \k\\ and |fc2| - * oo. There is only one exception where the integral converges: the segment with 6 = and 
7/2 < a < 4. This segment corresponds to wave action independent of vertical wavenumbers, dn/dm = 0. Within 
this segment we numerically determine a new steady convergent solution to Eq. ([1]) , with 

n(k,m) oc\k\- 3 - 7 . (3) 

This solution is not far from the large-wavenumber form of the Garrett-Munk (GM) spectrum: 

n(k,m) oc \k\~ 4 . (4) 

Alternatively, one can explore the physical interpretation of divergent solutions. We find a region in (a, b) space 
where there are both IR and UV divergences having opposite signs. This suggests a possible scenario where the two 
divergent contributions may cancel each other, yielding a steady state. An example of such a case is provided by the 
Pelinovsky-Raevsky (PR) spectrum, 

n k , m oc |fc|- 7 /a| m |-Va. (5) 

This solution, however, is only one among infinitely many. The problem at hand is a generalization of the concept 
of principal value integrals: for a and b which give opposite signs of the divergences at zero and infinity, one can 
regularize the integral by cutting out small neighborhoods of the two singularities in such a way that the divergences 
cancel each other and the remaining contributions are small. Hence all the exponents which yield opposite-signed 
divergence at the both ends can be steady solutions of Eq. |T]). As we will see below this general statement helps to 
describe the experimental oceanographic data which are available to us. The nature of such steady solutions depends 
on the particular truncation of the divergent integrals. 

So far, we have kept the formalism at the level of the self-similar limit of the kinetic equation |T]). However, once one 
considers energy transfer mechanisms dominated by interactions with extreme modes of the system, one can no longer 
neglect the deviations from self-similarity near the spectral boundaries: the inertial frequency due to the rotation of 
the Earth at the IR end, and the buoyancy frequency and/or dissipative cut-offs at the UV end. 

For example, we may consider a scenario in which interactions with the smallest horizontal wavenumbers dominate 
the energy transfer within the inertial subrange, either because the collision integral at infinity converges or because 
the system is more heavily truncated at the large wavenumbers by wave breaking or dissipation. We will demonstrate 
that the IR divergence of the collision integral has a simple physical interpretation: the evolution of each wave is 
dominated by the interaction with its nearest neighboring vertical wavenumbers, mediated by the smallest horizontal 
wavenumbers of the system, a mechanism denoted Induced Diffusion in the oceanographic literature. 

To bring back the effects of the rotation of the Earth in Eq. |T]), one introduces the Coriolis parameter / there and 
in the linear dispersion relation. Since we are considering the evolution of waves with frequency u much larger than 
/, / can be considered to be small. However, since the interaction with waves near / dominates the energy transfer, 
one needs to invert the order in which the limits are taken, postponing making / zero to the end. This gives rise to an 
integral that diverges like / raised to a negative power smaller than —1, but multiplied by a prefactor that vanishes 
if either 9 — 2a — 36 = or b = 0. These are the induced diffusion lines of steady state solutions, found originally in 
Ref. [H. This family of stationary states does a reasonable job of explaining the gamut of observed variability. 

The paper is organized as follows. Wave turbulence theory for the internal wave field and the corresponding kinetic 
equation are briefly summarized in Sec. [Til along with the motivating observations. We analyze the divergence of the 
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kinetic equation in Sec. IIIIl Section IIVI includes a special, convergent power-law solution that may account for the 
GM spectrum. In Sec. [V] we introduce possible quasi-steady solutions of the kinetic equation which are based on 
cancellations of two singularities. Section I VII shows that the IR divergence is dominated by induced diffusion, and 
computes the family of power-law solutions which arises from taking it into account. We conclude in Sec. IVIll 

II. WAVE TURBULENCE THEORY FOR INTERNAL WAVES 
A. Background and history 

The idea of using wave turbulence formalism to describe internal waves is certainly not new; it dates to Kenyon, 
with calculations of the kinetic equations for oceanic spectra presented in Refs. l6l lTolfl2l . 

Various formulations have been developed for characterizing wave- wave interactions in the stratified wave turbulence 
in the last four decades (see Ref. [13 for a brief review, and Refs. [l 2. 3. 4 5.6. 7.8.9. 10 illfl2lfl3l for details). 

We briefly discuss the derivation of the kinetic equation and wave-wave interaction matrix elements below in 
Eq. (HU). 

The starting point for the most extensive investigations has been a non-canonical Hamiltonian formulation in 
Lagrangian coordinates^ that requires an unconstrained approximation in smallness of wave amplitude in addition to 
the assumption that nonlinear transfers take place on much longer time scales than the underlying linear dynamics. 
Other work has as its basis a formulation in Clebsh-like variables^ and a non-Hamiltonian formulation in Eulerian 
coordinates 1 . Here we employ a canonical Hamiltonian representation in isopycnal coordinates^ which, as a canonical 
representation, preserves the original symmetries and hence conservation properties of the original equations of motion. 

Energy transfers in the kinetic equation are characterized by three simple mechanisms identified by Ref. and 
reviewed by Ref. Il8l. These mechanisms represent extreme scale-separated limits. One of these mechanisms represents 
the interaction of two small vertical scale, high frequency waves with a large vertical scale, near-inertial (frequency 
near /) wave and has received the name Induced Diffusion (ID). The ID mechanism exhibits a family of stationary 
states, i.e. a family of solutions to Eq. @. A comprehensive inertial-range theory with constant downscale transfer 
of energy can be obtained by patching these mechanisms together in a solution that closely mimics the empirical 
universal spectrum (GM). A fundamental caveat from this work is that the interaction time scales of high frequency 
waves are sufficiently small at small spatial scales as to violate the assumption of weak nonlinearity. 

In parallel work, Pelinovsky and Raevsky— derived a kinetic equation for oceanic internal waves. They also have 
found the statistically steady state spectrum of internal waves, Eq. ([5]), which we propose to call Pelinovsky- Raevsky 
spectrum. This spectrum was latter found in Refs. [l|5|. It follows from applying the Zakharov-Kuznetsov conformal 
transformation, which effectively establishes a map between the neighborhoods of zero and infinity. Making these two 
contributions cancel point- wisely yields the solution (O- 

Both Pelinovsky and Raevsky— and Caillol and Zeitlin— noted that the solution §5§ comes through a cancellation 
between oppositely signed divergent contributions in their respective collision integrals. A fundamental caveat is that 
one can not use conformal mapping for divergent integrals. Therefore, the existence of such a solution is fortuitous. 

Here we demonstrate that our canonical Hamiltonian structure admits to a similar characterization: power-law 
solutions of the form ^ return collision integrals that are, in general, divergent. Regularization of the integral allows 
us to examine the conditions under which it is possible to rigorously determine the power-law exponents (a, b) in 
Eq. ^ that lead to stationary states. In doing so we obtain the ID family. 

The situation is somewhat peculiar: We have assumed weak nonlinearity to derive the kinetic equation. The kinetic 
equation then predicts that nonlocal, strongly scale-separated interactions dominate the dynamics. These interactions 
have a less chance to be weakly nonlinear than regular, "local" interactions. Thus the derivation of the kinetic equation 
and its self-consistency is at risk. As we will see below, despite this caveat, the weakly nonlinear theory is consistent 
with much of the observational evidence. 



B. Experimental motivation 

Power laws provide a simple and intuitive physical description of complicated wave fields. Therefore we assumed that 
the spectral energy density can be represented as Eq. ([2]) , and undertook a systematic study of published observational 
programs. In doing so we were fitting the experimental data available to us by power-law spectra. We were presuming 
that the power laws offer a good fit of the data. We were not assuming that spectra are given by Garrett and Munk 
spectrum. That effort is reported in detail on elsewhere^, here we just give a brief synopsis to motivate present 
theoretical study. 

We reviewed the following observational programs: 
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FIG. 1: The observational points. The filled circles represent the Pelinovsky-Raevsky (PR) spectrum, the convergent numerical 
solution determined in Sec. IIVI (C) and the GM spectrum. Circles with stars represent estimates based upon one-dimensional 
spectra from the western North Atlantic south of the Gulf Stream (IWEX, FASINEX and SFTRE/PMIII), the eastern North 
Pacific (STREX/OS and PATCHEX 1 ), the western North Atlantic north of the Gulf Stream (Site-D), the Arctic (AlWEX) 
and the eastern North Atlantic (NATRE 1 and NATRE 2 ). There are two estimates obtained from two-dimensional data sets 
from the eastern North Pacific (SWAPP and PATCHEX 2 ) represented as circles with cross hairs. NATRE 1 represents a fit 
over frequencies of 1 < uj < 6cpd and NATRE 2 a fit over higher frequencies. 



• The Internal Wave Experiment (IWEX)2£; 

• The Arctic Internal Wave Experiment (AIWEX)2ii22; 

• The Frontal Air-Sea Interaction Experiment (FASINEX)22i2i; 

• The Patches Experiment (PATCHEX)^; 

• The Surface Wave Process Program (SWAPP) experiment; 

• North Atlantic Tracer Release Experiment (NATRE) 27 / SubductiorJS; 

• Salt Finger Tracer Release Experiment (SFTRE)^ / Polymode IIIc (PMIII) 30 ; 

• Site-D2i; 

• Storm Transfer and Response Experiment (STREX)^ / Ocean Storms Experiment (OS)^. 

Here we extract the summary of that study in the form of Fig. [TJ in which estimates of power-law exponents (a, b) 
are plotted against each other. 

That study comes with many caveats. The most important here is that the power-law exponents, (a, b), are 
not typically derived from two-dimensional horizontal wavenumber, vertical wavenumber spectra. The power laws 
were fit to one-dimensional frequency and one-dimensional vertical wavenumber spectra of energy, e(u>) oc lo~ c and 
e(m) ex; m~ d , with conversion to (a, b) requiring that: 

a = c + 2 and b — d — c. 



This procedure explicitly assumes that the two-dimensional spectra are separable. That is to say that the oj-m 
spectrum can be represented as a product of function of lo and function of m alone. It is clear that this assumption 
is not warranted in all cases. This is particularly true of the NATRE data, for which two power-law estimates are 
provided. Nevertheless, this is the best data that is available to us at a present time. We see that these points are 
not randomly scattered, but have some pattern. Explaining the location of the experimental points and making sense 
out of this pattern is the main physical motivation for this study. 
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C. Hamiltonian Structure and Wave Turbulence Theory 

This subsection briefly summarizes the derivation in Ref. [5|; it is included here only for completeness and to allow 
references from the core of the paper. 

The equations of motion satisfied by an incompressible stratified rotating flow in hydrostatic balance under the 
Boussinesq approximation are: 

d dz „ f dz \ 

V ■ — u =0, 



dt dp \ dp 

— + fu^ +u-Vu+ = 0, 

at po 

These equations result from mass convervation, horizontal momentum convervation and hydrostatic balance. The 
equations are written in isopycnal coordinates with the density p replacing the height z in its role as independent 
vertical variable. Here u — (u,v) is the horizontal component of the velocity field, u 1 - = (—v,u), V = (d/dx,d/dy) 
is the gradient operator along isopycnals, M is the Montgomery potential 

M = P + gpz, 

f is the Coriolis parameter, and po is a reference density in its role as inertia, considered constant under the Boussinesq 
approximation. 

The potential vorticity is given by 

/ + dv/dy — du/dx ^ 

where IT = p/gd 2 M/dp 2 — pdz/dp is a normalized differential layer thickness. Since both the potential vorticity and 
the fluid density are conserved along particle trajectories, an initial profile of the potential vorticity that is a function 
of the density will be preserved by the flow. Hence it is self-consistent to assume that 

where no(p) = —g/N(p) 2 is a reference stratification profile with the constant background buoyancy frequency, 
N = (-g/ipdz/dplbg)) 1 / 2 . This assumpt ion is not unrealistic: it represents a pancake-like distribution of potential 
vorticity, the result of its comparatively faster homogenization along than across isopycnal surfaces. 

It is shown in Ref. [5| that the primitive equations of motion @ under the assumption (jHJ) can be written as a pair 
of canonical Hamiltonian equations, 

sn sn d<p _ sn 

~dt ~ ' ~dt ~ ' ( ' 

where <f> is the isopycnal velocity potential, and the Hamiltonian is the sum of kinetic and potential energies, 



H = J dxdp(~(TL + IL(x,p)) 



V(f>(x,p) + ^-V ± A- 1 Il(x,p) 



dp -. 



Here, = (—d/dy,d/dx) and A 1 is the inverse Laplacian. 

Switching to Fourier space, and introducing a complex field variable a p through the transformation 



n p = n - 



¥1^1 
ivn |fc| 



V 2 9^v 



(a p + a*_ p ) 



(10) 



(11) 
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where the frequency uj satisfies the linear dispersion relation 



"v = \P + iU%, (12) 

V pfi" m 



the equations of motion ((6|) adopt the canonical form 



.d sn 

l m ap = ' (13) 



with Hamiltonian: 

H = J dpujp\a p \ 2 
+ I dp 012 (5 P 

+P1+P2 ,Pl,P2 a p a p 1 a p 2 "I C.C.J 

~r°-p+Pi+p 2 ,p 2 a p a Pi a P2 T C.C. 

)) • (14) 

This is the standard form of the Hamiltonian of a system dominated by three-wave interactions^. Calculations of 
interaction coefficients are tedious but straightforward task, completed in Ref. H These coefficients are given by 



N 1 

VpiP2 = ±jTgkk x k 2 ( T w» + J ^ pi + K P'Pi'Pi) ' 



Up ' PuP2 ~ 4^3 kkrk 2 ( Ip ' Pl ' P2 + Jp ? P2 + KpP1 ' 



(15a) 
(15b) 



r -- r lUJ2 k 2 k k 

1 p,pi,P2 — W ~ ™ "'l "'2 

- ((0, 1, 2) -> (1, 2, 0)) - ((0, 1, 2) -» (2, 0, 1)) , (15c) 
J p UP2 = - F f={k 2 k l -k 2 

- ((0, 1, 2) (1, 2, 0)) - ((0, 1, 2) - (2, 0, 1)) ) , (15d) 
Kp,p lt p 2 = — if\ ( \ [k\ — k 2 )k\ ■ k^ 



+ ((0,l,2)-f (1,2,0)) + ((0,1,2) -» (2, 0,1)) Jj, 



(15e) 



where ((0,1,2) — > (1,2,0)) and ((0,1,2) — > (2,0,1)) denote exchanges of suffixes. We stress that the field equation 
(|13|) with the three- wave Hamiltonian (1121 114| I15|) is equivalent to the primitive equations of motion for internal waves 
([6]). Other approaches, in particular the Lagrangian approach, is based on small-amplitude expansion to arrive to this 
type of equations. 

In wave turbulence theory, one proposes a perturbation expansion in the amplitude of the nonlinearity, yielding 
linear waves at the leading order. Wave amplitudes are modulated by the nonlinear interactions, and the modulation 
is statistically described by an kinetic equation^ for the wave action n p defined by 

n P S(p-p') = (a* p a p/ ). (16) 

The derivation of this kinetic equation is well studied and understood. For the three- wave Hamiltonian (|14p . the 
kinetic equation is the one in Eq. (TT]) , describing general internal waves interacting in both rotating and non-rotating 
environments. 

The delta functions in the kinetic equation ensures that spectral transfer happens on the resonant manifold, defined 
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P = Pl + P2 , (17a) 

Ld = LUl + L0 2 

(17b) 

P2=P + Pl . (17c) 

Now let us assume that the wave action is independent of the direction of the horizontal wavenumber, 

n p = n(\k\, \m\). 

Note that value of the interaction matrix element is independent of horizontal azimuth as it depends only on the 
magnitude of interacting wavenumbers. Therefore one can integrate the kinetic equation {T]) over horizontal azimutbM, 
yielding 

-Tjj- = I J (R12 - R lo - R Qi) dk 1 dk 2 dm 1 dm 2 , 
Rl2 = fpip 2 \Vpip2\ ^m-mi-TOa^Wp-Wpj -w P2 kk\k 2 / Si 2 ■ 

(18) 

Here S± 2 appears as the result of integration of the horizontal-momentum conservative delta function over all possible 
orientations and is equal to the area of the triangle with sides with the length of the horizontal wavenumbers k = \k\, 
k\ = I fci| and k 2 = \k 2 \. This is the form of the kinetic equation which will be used to find scale-invariant solutions 
in the next section. 

III. SCALE-INVARIANT KINETIC EQUATION 
A. Reduction of Kinetic Equation to the Resonant Manifold 

In the high-frequency limit to 3> /, one could conceivably neglect the effects of the rotation of the Earth. The 
dispersion relation (fT!?|) then becomes 

g \k\ 

cj p = u; fe , ro ~ — — r—r , (19) 
PoN \m\ 

and, to the leading order, the matrix element (| 15[) retains only its first term, ip, Pl , P2 . 

The azimuthally-integrated kinetic equation (|18|) includes integration over k\ and k 2 since the integrations over mi 
and m 2 can be done by using delta functions. To use delta functions, we need to perform what is called reduction to 
the resonant manifold. Consider, for example, resonances of type (|17a|) . Given k, k±, k 2 and m, one can find mi and 
m 2 satisfying the resonant condition by solving simultaneous equations 

k k\ k 2 

m = mi+m 2 , - — r = , r + r . (20) 

\m\ I mi I |m — mi| 

The solutions of this quadratic equation are given by 

mi = H (k + fci + k 2 + y/(k + h+ k 2 ) 2 - 4kk^ 

m 2 = m — mi., (21a) 
and 

mi = ^ (k - ki - k 2 - yj(k-ki - k 2 ) 2 + Akk^j 

m 2 = m - mi.. (21b) 
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Note that Eq. (|21ap translates into Eq. (|21b)) if the indices 1 and 2 are exchanged. Similarly, resonances of type (|17bl) 
yield 



(22a) 



m 2 


= -H(k 

2k V 


mi 


= m + m 2 


and 




m 2 


- -—(k 

- 2k r 


77li 


= m + m 2 



(k + ki - k 2 + ^Jik + kx -k 2 ) 2 +4kk 2 ^ 



(22b) 



and resonances of type (|17cp yield 



mi = 


2fc V 


m 2 = 


m + mi 


and 




mi = 


2fc V 


m 2 = 


m + mi 



(23a) 



(23b) 

After this reduction, a double integral over k\ and k 2 is left. The domain of integration is further restricted by the 
triangle inequalities 

k<ki + k 2 , ki<k + k 2 , and k 2 <k + ki. (24) 

These conditions ensure that one can construct a triangle out of the wavenumbers with lengths k, k\ and k 2 and 
determine the domain in the (fci, k 2 ) plane called the kinematic box in the oceanographic literature. 

Numerical evaluation of the collision integral is a complicated yet straightforward task. Interpretation of the results, 
though, is more difficult, mostly due to the complexity of the interaction matrix element and the nontrivial nature 
of the resonant set. Starting with McComas and Bretherton— , therefore, predictions were made based on a further 
simplification. This simplification is based on the assertion that it is interactions between wavenumbers with extreme 
scale separation that contribute mostly to the nonlinear dynamics. Three main classes of such resonant triads appear, 
characterized by extreme scale separation. These three main classes are 

• the vertical backscattering of a high-frequency wave by a low-frequency wave of twice the vertical wavenumber 
into a second high-frequency wave of oppositely signed vertical wavenumber. This type of scattering, as in 
Eqs. (|25al I27bl l29al [30a|) below, is called elastic scattering (ES). 

• The scattering of a high-frequency wave by a low-frequency, small-wavenumber wave into a second, nearly 
identical, high-frequency large- wavenumber wave. This type of scattering, as in Eqs. (125bl I27a1 129bl I30bj) below, 
is called induced diffusion (ID). 

• The decay of a small-wavenumber wave into two large vertical-wavenumber waves of approximately one-half 
the frequency. This type of scattering, as in Eqs. (|26al I26bl I28al I28b[) below, is called parametric subharmonic 
instability (PSI). 

To see how this classification appears analytically, we perform the limit of k\ — > and the limit fci — > oo in Eqs (|21l - 
|2"3"|) . We will refer to the k\ or k 2 — * limits as IR limits, while the fci and k 2 — * oo limit will be referred as an UV 
limit. Since the integrals in the kinetic equation for power-law solutions will be dominated by the scale-separated 
interaction, this will help us analyze possible solutions to the kinetic equation. 

The results of the k\ — > limit of Eqs. f2"lTl2"3")) are given by 

mi — > 2m, m 2 — > — m s 
-mi < m, m 2 — m fo^hl 

Wl«W,W 2 ~W 
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mi < to, m 2 ~ -m . , 

o j (26a) 

o i (26b) 



-mi < to, m 2 ~ m , . 

toi^-2to,to 2 ^-to , . 



We now see that the interactions (|25a[ 127b"]) correspond to the elastic scattering (ES) mechanism, the interactions (|25bl 
I27al correspond to the induced diffusion (ID). The interactions (|26al I26b|) . correspond to the parametric subharmonic 
instability (PSI). 

Similarly, taking the k\ and k 2 — > oo limits, of Eqs. (|2IH23p we obtain 



TOl 3> TO, — TO2 S> to 
— TOl 3> TO, TO 2 3> TO 



TOi ~ m/2, TO2 ~ —to/2 
wi, w 2 > w 

—mi 3> m, — to 2 ^> m 

mi ~ —m/2, to 2 ~ to/2 
— TOi 3> to, — m 2 ^> TO 



(28a) 
(28b) 

(29a) 
(29b) 

(30a) 
(30b) 



We now can identify the interactions (|28a[ I28b[) to be PSI, the interactions (|29al I30a|) to be ES, and finally the 
interactions (|29bl I30b[) as being ID. 

This classification provides an easy and intuitive tool for describing extremely scale-separated interactions. We will 
see below that one of these interactions, namely ID, explain reasonably well the experimental data that is available 
to us. 



B. Convergences and divergences of the kinetic equation. 



Neglecting the effects of the rotation of the Earth yields a scale-invariant system with dispersion relation given 
by Eq. (fT9)) and matrix element given only by the Ip, Pl . P2 in Eq. (|15p . This is the kinetic equation of Refs. I5ll3~il. 
describing internal waves in hydrostatic balance in a non-rotating environment. 

Proposing a self-similar separable spectrum of the form it is clear from the bi- homogeneous nature of the 
azimuthally-integrated kinetic equation (|18p that 

dn(ak,(3m) 4 -2a„i-2fe dn(k,m) 

dt = " P (31) 

In order to find a steady scale-invariant solution for all values of k and m, it is therefore sufficient to find exponents 
that give zero collision integral for one wavenumber. One can fix k and to, adopting for instance k — m = 1, and seek 
zeros of the collision integral as a function of a and b: 



(32) 
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FIG. 2: Divergence/convergence due to IR wavenumbers (top) and due to UV wavenumbers (bottom). The integral converge 
for the exponents in the shaded regions or on the segments (6 = 0,7/2 < a < 4) and (b = 1, 3 < a < 7/2). Dashed lines 
distinguish the domains where the indicated named triads dominate the singularity. 



Before embarking on numerical or analytical integration of the kinetic equation (|18| with scale-invariant solu- 
tions (p|. it is necessary to check whether or not the collision integral converges. Appendix lAl outlines these calcula- 
tions. The condition for the scale-invariant collision integral (|32[) to converge at the IR end, k\ or &2 — ► is given by 



or 



a + b/2- 7/2 < and - 3 < 6 < 3, 
a - 4 < and b = 0, 
a - 7/2 < and b= 1, 
a + b - 5 < and b > 3, 

a — 5 < and b < —3. 



(33a) 
(33b) 
(33c) 
(33d) 

(33e) 



Similarly, UV convergence as k\ and &2 - > oo implies that 

a + 6/2 - 4 > and - 2 < b < 2, (34a) 

a- 7/2 > and 6 = 0, (34b) 

a - 3 > and 6 > 2, (34c) 

or 

a + b - 3 > and 6 < -2. (34d) 
The domains of divergence and convergence are shown in Fig. [51 
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FIG. 3: Value of the collision integral as a function of a on the convergent segment, 6 = 0. 

Figure [2] also displays the classes of triads dominating the interactions. Knowing the classes of interactions that 
lead to the divergences of the kinetic equation allows us to find possible physical scenarios of the convergent solutions 
or to find a possible physical regularization of the divergences. 

Note that in addition to the two-dimensional domain of IR convergence [the regions (|33al I33dl I33c|) ] there are two 
additional IR convergent line segments given by Eqs. (|33b[) and (|33cp . These two special line segments appear because 
of the 6(6 — 1) prefactor to the divergent contributions to the collision integral (1A3|) . Similarly, for the UV limit, in 
addition to the two-dimensional region of convergence (Mil EU ESI) there is an additional special line segment of 

6 = piL) . 

We see that these domains of convergence overlap only on the segment 

7/2 < a < 4 and 6 = 0. (35) 

Note that 6 = corresponds to wave action independent of vertical wavenumbers, dn/dm = 0. Existence of the 6 = 
line will allow us to find novel convergent solution in Sec. HVl We also note that the IR segment on 6 = coincides 
with one of the ID solution determined in Sec. IVIl The other segment on 6 = 1 do not coincides with the ID solution 
in Sec. IVIl since the scale-invariant system has higher symmetry than the system with Coriolis effect. 

IV. A NOVEL CONVERGENT SOLUTION 

To find out whether there is a steady solution of the kinetic equation along the convergent segment ([33]) . we 
substitute the power-law ansatz ([2]) with 6 = into the azimuthally- integrated kinetic equation (fl"8|) . We then 
compute numerically the collision integral as a function of a for 6 = 0. To this end, we fix k = m = 1 and perform 
a numerical integration over the kinematic box (|24p . reducing the integral to the resonant manifold as described in 

sec.nng 

The result of this numerical integration is shown in Fig. [3] The figure clearly shows the existence of a steady 
solution of the kinetic equation (fT5|) near a = 3.7 and 6 = 0. 

This is, therefore, the only convergent steady solution to the scale-invariant kinetic equation for the internal wave 
field. It is highly suggestive that it should exist so close to the GM spectrum, a — 4 and 6 = for large wavenumbers, 
the most agreed-upon fit to the spectra observed throughout the ocean. It remains to be seen whether how this 
solution is modified by inclusion of background rotation. 

V. BALANCE BETWEEN DIVERGENCES 

The fact that the collision integral C diverges for almost all values of a and 6 can be viewed both as a challenge 
and as a blessing. On the one hand, it makes the prediction of steady spectral slopes more difficult, since it now 
depends on the details of the truncation of the domain of the integration. Fortunately, it provides a powerful tool 
for quantifying the effects of fundamental players in ocean dynamics, most of which live on the fringes of the inertial 
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FIG. 4: Signs of the divergences where the both IR and UV contributions diverge. The left symbols show the signs due to IR 
wavenumbers and the right symbols show the signs due to UV wavenumbers. 



subrange of the internal wave field: the Coriolis effect, as well as tides and storms at the IR end of the spectrum, and 
wave breaking and dissipation at the UV end. The sensitive response of the inertial subrange to the detailed modeling 
of these scale-separated mechanisms permits, in principle, building simple models in which these are the only players, 
bypassing the need to consider the long range of wave scales in between. 

At the IR end, the resonant interactions are dominated by the ID singularity for —3 < b < 3, Fig. O The sign of 
the divergences is given by —b(b— 1) (Eq. (|A3[) V Similarly, at the UV end resonant interactions are dominated by the 
ID singularity for —2 < b < 2 (Fig. [2]). The sign of the divergence is given by b (Eq. (|A4[) ). At the UV end for b > 2, 
where ES determines the divergences, the sign of the singularity is given by —b, i.e. the sign is negative (Eq. (|A5|) ). 
Figure 2] shows the signs of the divergences where both the IR and UV contributions diverge: the left sign corresponds 
to the IR contribution, and the right sign to the UV contribution. 

Hence, in the regions, 

7<2a + 6<8, and - 2 < b < 1, (36a) 

or 

7 < 2a + 6, a < 3, and b > 2, (36b) 

the divergences of the collision integral at the IR and UV ends have opposite signs. Then formal solutions can be 
constructed by having these two divergences cancel each other out. 

This observation justifies the existence of the PR solution ([5]). Indeed, the PR spectrum has divergent power- 
law exponents at the both ends. One can prove that the PR spectrum is an exact steady solution of the kinetic 
equation |T]) by applying the Zakharov-Kuznetsov conformal mapping for systems with cylindrical symmetr y 35 i 36 ' 37 . 
This Zakharov-Kolmogorov conformal mapping effectively establishes a map between the neighborhoods of zero and 
infinity. Making these two contributions cancel point-wisely yields the solution ([5]). For this transformation to be 
mathematically applicable the integrals have to converge. This transformation leads only to a formal solution for 
divergent integrals. 

The PR spectrum was first found by Pelinovsky and Raevsky 11 . However, they realized that it was only a formal 
solution. The solution was found again in Ref.[l] through a renormalization argument, and in Ref.|H within an isopycnal 
formulation of the wave field. 

The idea of a formal solution, such as PR, can be generalized quite widely: in fact, any point in the regions with 
opposite-signed divergences can be made into a steady solution under a suitable mapping that makes the divergences 
at zero and infinity cancel each other, as does the Kuznetsov-Zakharov transformation. Such generalized Kuznetsov- 
Zakharov transformation is an extension of the idea of principle value for a divergent integral, whereby two divergent 
contributions are made to cancel each other through a specific relation between their respective contributions. 

This scenario also suggests an possible explanation for the variability of the power-law exponents of the quasi-steady 
spectra: inertial subrange spectral variability is to be expected when it is driven by the nonlocal interactions. The 
natural local variability of players outside the inertial range, such as tides and storms, translates into a certain degree 
of non-universality. 
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VI. REGULARIZATION BY THE CORIOLIS EFFECT 



Physically, ocean does not perform generalized Kuznetsov-Zakharov transformation. However, in the ocean there 
are finite boundaries in frequency domain. In particular, the inertial frequency, /, provides a truncation for the IR 
part of the spectrum, while the UV truncation is provided by the buoyancy frequency, N. These two frequencies vary 
from place to place, giving grounds for spectral variability. Consequently, the integrals are not truly divergent, rather 
they have a large numerical value dictated by the location of the IR cutoff. 

All the experimental points are located in regions of the (a, b) domain a > (— b + 7)/2 for which the integral diverges 
in the IR region. For most of the experimental points, the UV divergence is not an issue, as a > —6/2 + 4. The UV 
region is therefore assumed to be either sub-dominant or convergent in this section, where we study the regularization 
resulting from a finite value of /. 

Since the IR cutoff is given by /, a frequency, it is easier to analyze the resulting integral in (cji, mi) rather than in 
the traditional (&i,mi) domain. Thus we need to express both the kinetic equation and the kinematic box in terms 
of frequency and vertical wavenumber. For this, we use the dispersion relation (|12|) to express k in terms of lo in the 
description of the kinematic box 



u>i < E 3 (mi), uji > E 4 (mi), 


loi > Ei(mi) 




if mi < 0, u)\ > lo, 




(37a) 


oji > E 3 (mi), uji < S 4 (mi), 


loi < Ei(mi) 




if mi < 0, lo\ < lo, 




(37b) 


u>i < £2 (mi), loi > E 2 (mi), 


loi > £4 (mi) 




if mi > 0, loi < lo, 




(37c) 


loi > E a (mi), loi < Ei(mi), 


loi < E 2 (mi) 




if mi > 0, loi > lo, 




(37d) 



where we have introduced the four curves in the (cji,mi) domain that parameterize the kinematic box: 

Ei(loi) = 



-V-f 2 


+ L0 2 


+ V-f 2 


+ {lo 


-Ul) 2 


v-p- 


hLof H 


-V-f 2 - 


- (lo - 


~^l) 2 


v 7 -/ 2 + 


LO 2 + 


V-f 2 + 


(lo- 


LOl) 2 


V-f 2 + 


LO 2 + 


V-f 2 + 


(lo- 


LOl) 2 ' 


-V-f 2 


+ L0 2 


+ V-f 2 


+ [lo 


-001) 2 


-V-f 2 


+ L0'f 


+ V-f 2 


+ {lo 


-LOl) 2 


V-.f 2 


\- LO 2 H 


- V-f 2 H 


- (lo - 


-LOl) 2 


-V-f 2 


+ LO~( 


+ V-P 


+ (lo 


-LOl) 2 



E 2 (loi) = 
E 3 (loi) = 
E 4 (loi) = 

The kinematic box in the (lo, m) domain is shown in Fig. [5] To help in the transition from the traditional kinematic 
box to the kinematic box in (LO,m) domain, the following limits were identified: 

• ID1 is the ID limit of Eq. (|25b[) with indices 1 and 2 being flipped 

• ID2 is the ID limit of Eq. (|27a|) with indices 1 and 2 being flipped 

• ID3 is the ID limit of Eq. (|25bl 

• ID4 is the ID limit of Eqs. (f29bl l30bj) 

• PSI1 is the PSI limit of Eq. (|28a|) 

• PSI2 is the PSI limit of Eq. ((28b)) 

• PSI3 is the PSI limit of Eq. ([26a|) 

• PSI4 is the PSI limit of Eq. ([26b]) 

• ESI is the ES limit of Eq. (|25a| 

• ES2 is the ES limit of Eq. ()27b[) with indices 1 and 2 being flipped 



14 



W-f Ul + f 



3 
2 



-1 
-2 



i PSI1 






Ei 

e\ 




- 


" ES1 

ID1, 


ID2 






- 
- 






PS 3 




ESI 






PSI4 ._ 








,/ES2 








ES3 


















ID4 


PSI2 











12 3 4 



FIG. 5: The kinematic box in the (wi,mi) domain. Two disconnected regions where uji < uj depict regions with "sum" 
interactions, namely uj = u)i + 01%, m = mi + rri2, type of the resonances. The connected regions where wi > u) depict 
"difference" resonances, cu2 = uJi — UJ„ mi = mi — m. The parameters are chosen so that f /uj = 0.1. 



• ES3 is the ES limit of Eq. (|27bj) with indices 1 and 2 being flipped 

• ES4 is the ES limit of Eq. (j29a|) 

An advantage of the {uj, m) presentation for the kinematic box is that it allows a transparent reduction to the resonant 
manifold. A disadvantage is the curvilinear boundaries of the box, requiring more sophisticated analytical treatment. 
To proceed, we assume a power-law spectrum, similar to Eq. ([2]), but in the {uj, m space: 

n uj,m °f uj a m b . (38) 

We need to transform the wave action as a function of k and m to a function of uj and m. This is done in Appendix IB"1 
The relation between a, b and a, b reads: 

a = —a, b = —a — b. 

Equation (JTSJ) transforms into 



d If \V° I 

—n{k{uj,m),m) = — I duj\dmiJ q 2 {nin 2 — n{n\ + n 2 )) 
Ov K J 12 

2 r iv 1 1 2 

— - / duJidmiJ 02 {nn 2 — ri\{n + n 2 )) 
k J b 2 



uj2—u)—uj\ ,m2—m—mi 



u)2—^\ -w,m2-mi-m 



(39) 



We have used the dispersion relation ki = rriiWuj'f — f 2 , and defined J as the Jacobian of the transformation from 
{k\,k 2 ) into {uji,uj 2 ), times the kk±k 2 factor, 



J = kk\k 2 



dk\ dk 2 



' dui du 2 

In Fig. there are three ID corners with significant contribution to the collision integral: 

1. ID1 region. In this region, mi is slightly bigger than m, uj\ slightly smaller than uj, and uj 2 and m 2 are both 
very small. This region can be obtained from the region l|37b[) above by interchanging indices 1 and 2. In this 
region, 

n 2 » n, n\. 

2. ID2 region. In this region, uj\ is slightly bigger than uj + f where uj 2 — uj± — uj is small, and m 2 = m\ — m is 
negative and small. This is the region (|37d[) . Also 



n 2 » n, m. 
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3. ID3 region. Small u>i, small negative mi. This corresponds to the region (|37bp . where 

n-i > n, ri2- 

Note that this region can be obtained from the region ID1 by flipping indices 1 and 2. Consequently, only one 
of the ID1, ID2 should be taken into account, with a factor of two multiplying the respective contribution. 

Making these simplification, and taking into account the areas of integration in the kinematic box, we obtain 

— n(k(uj,m),m) = - I dwi / draiJ— — ni{ni - n) - - I du>i I dmiJ—^ — m(n-ni), 

(40) 

where the small parameter to s is introduced to restrict the integration to a neighborhood of the ID corners. The 
arbitrariness of the small parameter will not affect the end result below. 

To quantify the contribution of near-inertial waves to a (w, ni) mode, we write 

e~ / <u; = 1. 

Subsequently, near the region ID3 of the kinematic box, we write 

Wi = / + e, 

while near the ID2 corners of the kinematic box we write 

u>i = cj + / + e. 

We then expand the resulting analytical expression (|40p in powers of e and / without making any assumptions on 
their relative size. These calculations, including the integration over vertical wavenumbers mi, are presented in 
Appendix [Cj The resulting expression for the kinetic equation is given by 

^n(k(uj,m),m) = ^- ~tj (5 - 3 (s + bj^ 

xto 5+2 V 3 + 5 -^ r de(e + /) 4+H+ '(e 2 + 2ef + 17 f) . 
Jo 

(41) 

The integral over e diverges at e = 0, if / = and if6 + a + 6>— 14^ 

However, if we postpone taking / = limit, we see that the integral is zero to leading order if 

2-3^3 + 6^=0 or a-6 = (42) 

or, in terms of a and 6, 

9 - 2a - 3b = or 6 = 0. (43) 

This is the family of power-law steady-state solutions to the kinetic equations dominated by infra-red ID interactions. 
These steady states are identical to the ID stationary states identified by McComas and Bretherton— , who derived 
a diffusive approximation to their collision integral in the infra-red ID limit. Note that McComas and Miiller— 
interpreted 6 = as a no action flux in vertical wavenumber domain, while 9 — 2a — 36 = is a constant action flux 
solution. What is presented in this section is a rigorous asymptotic derivation of this result. These ID solutions helps 
us to interpret observational data of Fig. Q] that is currently available to us. 

VII. CONCLUSIONS 

The results in this paper provide an interpretation of the variability in the observed spectral power laws. Combining 
Figs. [HE] and H with Eqs. (4}, © and (43]), produces the results shown in Fig [6] 
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horizontal power-law exponents of 3D action spectrum (a) 



FIG. 6: The observational points and the theories. The filled circles represent the Pelinovsky-Raevsky (PR) spectrum, the 
convergent numerical solution determined in Section 4 and the GM spectrum. Circles with stars represent power-law estimates 
based upon one-dimensional spectra. Circles with cross hairs represent estimates based upon two-dimensional data sets. See 
Fig. [Tjfor the identification of the field programs. Light grey shading represents regions of the power-law domain for which the 
collision integral converges in either the IR or UV limit. The dark grey shading represents the region of the power-law domain 
for which neither the IR or UV limits converge. The region of black shading represents the sub-domain for which the IR and 
UV divergences have the same sign. Overlain as solid white lines are the induced diffusion stationary states. 

Notice that most of the observational points are UV convergent and IR divergent and that regularization of the 
collision integral by /, as in Sec. IVI1 produces a family of stationary states that collapses much of the variability. 
Four points land in the region where the spectra are both IR and UV divergent, with all but one point close to the ID 
lines. None of the observations lies in a region where both IR and UV divergences have the same sign. Furthermore, 
the novel convergent (non-rotating) solution is in close proximity to the experimental points determined from the 
two-dimensional spectra. 

Summarizing the paper, we have analyzed the scale-invariant kinetic equation for internal gravity waves, and shown 
that its collision integral diverges for almost all spectral exponents. Figure [6] shows that the integral nearly always 
diverges, either at zero or at infinity. This means that, in the wave turbulence kinetic equation framework, the energy 
transfer is dominated by the scale-separated interactions with either large or small scales. 

The only exception where the integral converges is a segment of a line, 7/2 < a < 4, with 6 = 0. On this convergent 
segment, we found a special solution, (a, b) = (3.7,0). This new solution is not far from the large-wavenumber 
asymptotic form of the Garrett-Munk spectrum, (a,b) = (4,0). 

We have argued that there exist two regions of power-law exponents which can yield quasi-steady solutions of the 
kinetic equation. For these ranges of exponents, the contribution of the scale-separated interactions due to the IR 
and UV wavenumbers can be made to balance each other. The Pelinovsky-Raevsky spectrum is a special case of this 
scenario. None of the observational power-law exponents lie within the region of same-signed UV and IR divergences. 

This scenario, in which the energy spectrum in the inertial subrange is determined by the nonlocal interactions, 
provides an explanation for the variability of the power-law exponents of the observed spectra: they are a reflection 
of the variability of dominant players outside of the inertial range, such as the Coriolis effect, tides and storms, can 
not be universal. 

We also rederive ID steady-solution lines, which are also shown as white lines in Fig. [5] 

All of theory, experimental data, and the results of numerical simulation in Ref. |38| hint at the importance of the 
IR contribution to the collision integral. The nonlocal interactions with large scales will therefore play a dominant 
role in forming the internal-wave spectrum. To the degree that the large scales are location dependent and not 
universal, the high-frequency, high vertical-wavenumber internal- wave spectrum ought to be affected by this variability. 
Consequently, the internal-wave spectrum should be strongly dependent on the regional characteristics of the ocean. 
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TABLE I: Asymptotics as ki -» 0. ES (pTal [23b]) gives e" a+5 owing to the symmetry of y. ID (f21bl [23a)l gives e - a -( 6 - 7 )/ 2 
(e _a+4 ) owing to the second cancellation. PSI (|22al I22b|) gives e - a ~ b + 5 ^ The asymptotics for 6 = appear in parentheses. 
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APPENDIX A: ASYMPTOTICS OF COLLISION INTEGRAL IN INFRA-RED AND ULTRA-VIOLET 

LIMITS 



Let us integrate Eq. (|18|) over mi and mi- 

drip 
~~dt 

rpO _ 

g°[ 2 (k 1 ,k 2 ) = 



3i, 2 (wi) = — 



k 2 



where g®' 2 appears owing to <5, 



m \mi | \m — m\ | ' 
and m\{k\,k<x) is given by the resonant conditions ([2"Tj) - lj2"3")) . 



(Al) 



1. Infra-red asymptotics 

We consider the asymptotics of the integral in Eq. (|A1[) as k% — * 0. We employ the independent variables x and 
y, where k\ = kx, k 2 — k(l + y), x,y — 0(e), x > and —x < y < x. In this limit of e — > 0, ni S> n,n 2 . In this 
limit, Eqs. (j21a)) and (j23bf . Eqs. (f2lE|) and (l23a|) . and Eqs. ([22a| and ([22b]) correspond to ES, Eqs. (j25al [2Tb|) ID, 
Eqs. (|25bl 127a) . PSI, Eqs. (|25al I27b]l . respectively. Without loss of generality m is set to be positive. 

Assuming the power-law spectrum, n(fc, m) — |fe| _0 |m| _b , we make Taylor expansion for the integrand of the kinetic 
equation (|A1[) as powers of e, that is x and y. Then, we get Table fl] which shows the leading orders of the each terms 
according to the asymptotics. The leading order of the collision integral is given by ID when —3 < b < 3. Therefore, 
we are going to show the procedure to get the leading order of ID (|21b[) and (|23a[) below. 

As e — > 0, n 2 — > n for ID solutions. Therefore, the leading orders of /f 2 ~ ni(n 2 — n) and / 2 j ~ ni(n — n 2 ) is 
(9( e - a -( fc - 1 )/ 2 ^ The order 0(e~ a ~ b / 2 ) is canceled as e — ► 0. This is called the first cancellation. It must be noted 
that the leading order when b = is 1/2 larger than that when 6^0 since dn/dm = 0. The leading orders when 
6 = are written in parentheses in Table [T] 

The leading order of the integrand in Eq. (|A1[) is written as 

-a-(h+l)/2 

T?,2-Tl -Tl lO tk- 2a+3 m- 2b+1 V 



y/(x + y)(x-y) 
x (-2ay 2 - b ((1 - b)y(x + y) - 2x(x - y)) + b(b + l)xy) . 

(A2) 
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TABLE III: Asymptotics as fei -» oo. PSI pTa ^ l21b ]l gives e a+b - 3 . ES J22al [23aJ| gives e a ~ 3 . ID (f22bl [23b]) gives e a + 6 / 2 - 4 
(e a_7 ' 2 ). The asymptotics for b — appear in parentheses. 
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Therefore, the integrand has 0(e a ^ 5 ^ 2 ). The term which has 0(e a b / 2 + 2 ) is canceled since XJ 2 x — > T{* 2 (and 
T 2 1 — ► 0) as e — > 0. This is the second cancellation. 

Finally, we get the leading order of the kinetic equation after integration over y from — x to x: 

^ oc -b{l-b)k 4 - 2a m 1 - 2b J^x- a -^' 2 dx. (A3) 

The integral has 0(e~ a ~( b ~ 7 ^ 2 ). Consequently integral converges if 

a + (b - 7)/2 < and - 3 < b < 3. 

The integral for the PR spectrum, which gives 0(e -1 / 4 ), diverges as k\ — > 0. However, the integral for the GM 
spectrum converges because 6 = and the next order is 0(1). 

It should be noted that the leading order when b = 1 is 1/2 larger than that when b ^ 0, 1 since a balance between 
first- and second-order derivative is made. The leading orders when b = 1 arc 0(e~ a+7 / 2 ). It is also helpful to note 
that T° 2 — T^q — T 2 Q = 0(e~ a+2 ) for ES because of no second cancellation. However, the collision integral has 
0(e~ a+5 ) because of symmetry of y. Therefore, the integral which is dominated by ES converges 

a — 5 < and b < —3. 

Similarly, the integral which is dominated by PSI converges 

a + b - 5 < and b > 3. 



2. Ultra-violet asymptotics 

Next, we consider the limit k\ —>■ oo. In this case, ki also approaches to infinity. We employ the independent 
variables x and y as ki = fc/2(l + l/x + y) and fc 2 = fc/2(l + l/a; — y), where x = O(e) and — 1 < y < 1. Again, m > 
is assumed. 

The leading orders are obtained by the similar manner used in the IR asymptotic and are summarized in Table HTT1 
The leading order of the integral is given by ID, whose wavenumbers are given by Eqs. (|22bl 123b)) . when —2 < b < 2. 
In this limit, no second cancellation is made. 

As the result of the perturbation theory, we get the leading order, 

^ oc k^m^b J^x a+b ^- 5 dx. (A4) 

It has 0(e a+b / 2 - 4 ). Therefore, the integral converges if 

a + b/2 - 4 > and - 2 < b < 2. 

The integral for the PR spectrum, which gives 0(e-^ 4 ), diver ges as k\ — > oo. and that for the Garrett-Munk 
spectrum, which gives O(e ), converges owing that 6 = 0. 
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dUp oc -Jfe^m 1 - 2 ^ f x a ~ 2 dx (A5) 
Jo 



Similarly, 

for ES, which is dominant for b > 2. Consequently, the integral converges also if 

a - 3 > and 6 > 2. 
In the same manner, the convergent domain of the integral for PSI is given by 

a + b - 3 > and b < -2. 

APPENDIX B: FREQUENCY-VERTICAL- WAVENUMBER AND 
HORIZONTAL- VERTICAL- WAVENUMBER SPECTRUM 

The theoretical work presented below addresses the asymptotic power laws of a three-dimensional action spectrum. 
In order to connect with that work, note that a horizontally isotropic power-law form of the three-dimensional wave 
action n(k, m) is given by Eq. j2]). 

The corresponding vertical wavenumber-frequency spectrum of energy is obtained by transforming rifc jTO from 
wavenumber space (k,m) to the vertical wavenumber-frequency space (cj,m) and multiplying by frequency. In the 
high-frequency large-wavenumber limit, 

E(m,w) oc uj 2 - a m 2 - a - b . 
The total energy density of the wave field is then 

E = J oj(k,Tn)n(k,m) dkdm = J E(cj,m) dujdm. 

Thus, we also it convenient to work with the wave action spectrum expressed as a function of to and m Therefore we 
also introduced (f3"5| . The relation between a, b and a, b reads: 

a = —a, b = —a — b. 



APPENDIX C: ASYMPTOTIC EXPANSION FOR SMALL / VALUES. 



In this section we perform the small / calculations of Sec. IVII We start from the kinetic equation written as 
Eq. (TiO]) . There we change variables in the first line of Eq. (|40[) as 

mi = £3(011) + /u(£i(a>i) - £ 3 (wi)), 

and in the second line of Eq. (|40|) as 

mi = E 3 (u>i) + /Lt(£ 2 (wi) - £ 3 (wi)). 
Then the Eq. (TiD|) becomes the following form: 

—n(k(LU,m),m) = — / dui\ I d[iP\ 
at k If In 



U-f pi 

dui / dnV 2 - (CI) 

u — f — u s JO 



Here we introduced integrand V\ and Vi to be 



l^i° 2 | 2 

^1 = J—^ — ni(n 2 - n)(£i(wi) - £ 3 (wi)), 
"1,2 

\V° I 2 

^2 = J-^§-n 2 {n - ni)(£ 2 (o;i) - £3(01)). (C2) 
"1,2 
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Before proceeding, note the following symmetry: 

E\(u)\ — lj — u>i) — m — E2{uo'i), 

and 

Eq(ll>i = LJ — Lu[) — JTl — E^LHi). 

To quantify the contribution of near-inertial waves to a (w, m) mode, we write 

e~ / < w = 1. 

Subsequently, in the domain (a) we write 

oji = f + e, 

in Vx, and 

0J\ = uj + / + e 

in Vi. Furthermore, we expand V\ and V2 in powers of e and / without making any assumptions of the relative 
smallness of / and e. We use the facts that 

m > 0, e > 0, / > 0, < fi < 1. 

Define 

Vi = P i +P il 

and 

P 2 = Pz + Pi- 

This allows us to expand Pi, P2, P3 and P4 in powers of / and e. We perform these calculations on Mathematica 
using Series command, and extensively using Assumptions field in the FullSimplify command. 

Mathematica was then able to perform the integrals of V\ and V2 over \i from to 1 in (|C1|) analytically. The 
result is given by Eq. (|4Tj) . 
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